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ABSTRACT 


The coordination of multifunction phased array radars across networked platforms can en¬ 
able superior functionality and battle space awareness. This thesis formulates and solves 
a number of optimization models and heuristic algorithms to analyze and prescribe radar- 
to-target assignments and schedules. One optimization model uses full target information 
to provide a “best case” assessment of the ability of a given set of radar platforms to track 
a collection of targets. A modified version of this model determines the impact on these 
results if targets coordinate their maneuvers in order to overwhelm the radar system. We 
then consider the more realistic scenario in which the planner’s knowledge is imperfect 
and describe approaches for allocating radar assets to targets in that setting. The first such 
approach extends an existing two-dimensional geographic allocation method to three di¬ 
mensions. This allows for an allocation of the operating space to radar assets and can 
serve as a preallocation heuristic for more sophisticated assignment algorithms. Moreover, 
because the existing method does not account for transfers in tasks between geographical 
areas, this thesis models the additional workload involved in performing handoffs of targets 
between radars. 
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Executive Summary 


Advanced warships cost in excess of one billion dollars yet remain vulnerable to a single 
missile. The proliferation of relatively low cost and increasingly sophisticated unmanned 
aerial vehicles, missiles, and decoys provides our enemies the capacity to overwhelm a 
ship’s defensive capabilities. The probability of damage from any one threat may remain 
low, but the success of a single missile could incapacitate a high value asset. The coordina¬ 
tion of multifunction phased array radars across platforms is critical to obtaining superior 
air defense and battle space awareness during periods of high contact density. This thesis 
develops several models to analyze target assignments in a multiple-sensor network. 

One optimization model provides a “best case” assessment of the ability of a given set of 
radar platforms to track a collection of targets. We assume the set of radars have perfect 
knowledge of the targets, including the times that each target track begins and ends and 
the location of the target throughout this period. This allows the model to optimize the 
schedule of radar observations over the entire time horizon and provides a bound for the 
real capabilities of a set of radars. A modified version of this model determines the impact 
on these results if targets coordinate their maneuvers in an attempt to exceed the resource 
capacity of the sensor network. Results suggest that synchronized maneuvers may produce 
situations where sensor resource demands exceed network capabilities. This could result 
in reduced track quality leading to higher track loss rates. 

We also consider the more realistic scenario in which the planner’s knowledge is imper¬ 
fect and describe approaches for allocating sensor assets to targets based on geographic 
location. The first method extends an existing two-dimensional geographic allocation ap¬ 
proach to three dimensions. By limiting the amount of the operating space we assign to 
a specific sensor, we develop a heuristic approach for preassigning sensors to targets that 
reduces the computational demands of more sophisticated assignment algorithms. In the 
scenarios examined, the heuristic can perform approximately 40 percent of sensor-to-target 
assignments with limited impact on the “best case” assessment results. This could lead 
to significant reductions in the real-time computations and data transfer rates required by 
sensor-to-target assignment algorithms. 
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Finally, we investigate the impact of transferring target tracking requirements between sen¬ 
sors by developing a model that permits the inclusion of handoff penalties. In order to do 
so, we modify the previously mentioned geographic allocation model by discretizing the 
operating space and defining target paths through this space. To account for the additional 
resources required when transferring a target between sensors, we assess a handoff penalty 
for any sensor that gains part of a target’s path after it begins. For the scenario consid¬ 
ered, we find that significant modifications to the sensors’ regional assignments only occur 
if the resource demands for target handoffs account for more than 25 percent of the total 
demand. 



CHAPTER 1: 
INTRODUCTION 


Multifunction phased array radars (MFRs) are capable of performing various tasks in rapid 
succession. The performance of target search, detect, and track operations concurrently 
with missile guidance functions allow MFRs to deliver superior battle space awareness 
and air defense capabilities than previously provided by multiple, single-purpose radars. 
Equipped with a communication network capable of sharing target data, a battle group can 
leverage the radar resources available on its ships to extend threat detection ranges, enhance 
track data, and synchronize defensive missile capabilities across platforms. 

While a typical shipborne MFR can track over 100 targets independently (U.S. Navy, 
2013), this could prove inadequate in a dense target environment. Combined with friendly 
air assets, the expanded availability of unmanned aerial vehicles (UAVs) and decoys pro¬ 
vides the opportunity to overwhelm the tracking capacity of a single MFR with large salvos 
and swarm tactics. MFRs require a coordinated approach to maximize search capabilities 
while maintaining track quality. The multifunctional nature of modern radars complicates 
this process, and platform-level demands or limitations on MFRs may reduce the capacity 
available for coordinated search and tracking. This thesis focuses on the allocation of radar 
tracking requirements across multiple MFRs to reduce radar resources required for track 
maintenance. In this thesis, we define any radar contact as a target, including neutral and 
friendly air assets, to conform with the established academic lexicon. 

1.1 Motivation and Background 

Advanced warships cost in excess of one billion dollars yet remain vulnerable to a single 
missile. The proliferation of relatively low cost and increasingly sophisticated UAVs, mis¬ 
siles, and decoys provides the capacity to overwhelm a ship’s defensive capabilities. The 
probability of damage from any one threat may remain low, but the success of a single mis¬ 
sile could incapacitate a high value asset (HVA). Coordination of radar resources within a 
battle group can increase tracking performance during periods of high contact density. The 
successful integration of radar resource management across platforms relies on two critical 
pieces of shipboard infrastructure: MFRs and tactical data links (TDLs). 
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1.1.1 Phased Array Radar 

MFRs consist of several elements, or antennas, spatially separated in an array. The array 
is designed to inerease antenna gain, to allow for eleetronie beam formation and steering, 
and to determine the direetion from whieh ineoming signals arrive. In passive MFRs, a 
single power souree provides energy to phase shifters attaehed to eaeh element. When the 
phase of emitted radiation between individual elements is shifted, eonstruetive interferenee 
patterns inerease energy propagated in the desired beam direetion while destruetive inter¬ 
ferenee reduees radiated energy emitted in other direetion. By altering the amount of phase 
shift between elements, the array ean direet radar beams eleetronie ally without the need to 
physieally move any eomponents, as illustrated in Figure 1.1. 



Figure 1.1: Beam steering with a linear phased array (from U.S. Navy Naval Air Systems 
Command, 2012, p. 3-4.3). 


Unlike traditional rotating dish radars, eleetronie beam steering allows an MFR to eyele 
through tasks in various direetions without delay. Thus, rather than performing “traek 
while sean” (TWS), where radar returns obtained during a seareh are filtered and eorrelated 
to form traeks, MFRs more eommonly operate in a “track and scan” (TAS) mode. In TAS, 
MFRs employ traeking beams, distinet from those generated for seareh, for eaeh target 
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being tracked. This allows the MFR to use advanced, higher energy techniques to improve 
the position accuracy of targets. TAS also allows the MFR to continue to maintain tracks 
outside the area in which it may be performing search functions. Additionally, TAS allows 
the MFR to maintain different data update rates for tracking and searching. For example, 
where a search area may need to be covered every 10 seconds, a known target may have to 
be observed every second (Billiter, 1989, pp. 89-92). 

While TAS allows for greater accuracy of target data, it also requires significantly more 
time and power from the radar than TWS. As the number of tracks increases, TAS may 
exceed the resources available to the MFR. The MFR can reduce this burden by using TWS 
to monitor long-range or low-priority targets while maintaining high-risk, high-priority 
targets in TAS (Billiter, 1989, pp. 93-96). 

To prioritize and perform its assigned tasks, an MFR utilizes a radar task scheduler (RTS). 
For each time interval Tp, the RTS accepts the following inputs: a set of search beams 
necessary to explore the assigned search region, a track task list with the required renewal 
time for each track, an auxiliary task list with renewal times for each task, and the time to 
execute each search, track, and auxiliary task. From this input, the RTS develops a task 
sequence for the MFR to perform over the time interval Tp and reports any modifications to 
input parameters it made while developing the task sequence (Sabatini & Tarantino, 1994, 
pp. 219-255). 

The electronic steering capability of MFRs permits the RTS to schedule several track tasks 
concurrently. For targets at an estimated position, the RTS can calculate the expected time 
between pulse transmission and a radar return. Rather than allow the MFR to idle during 
this delay, the RTS schedules additional tracking tasks such that the transmissions and 
returns for each target do not coincide. This process, known as pulse packing (Bogler, 
1990, pp. 252-253) or track interleaving (Billiter, 1989, pp. 99-101), allows an MFR to 
update several targets during the same time required to accomplish a single tracking task. 

MFRs have been in use for over three decades for aircraft surveillance and tracking (Na¬ 
tional Research Council, 2008, p. 19). Various models and variations are now installed on 
warships around the world. Most common on U.S. warships is the AN/SPY-1 radar. The 
AN/SPY-1 consists of four arrays positioned to allow 360 degree coverage. It transmits in 
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the 2 to 4 Ghz range with a peak power of 4 to 6 MW. As the primary sensor for the Aegis 
Weapon System, it eonducts horizon and sector searches concurrently with multiple target 
tracking and designation (Jane’s, 2012). Figure 1.2 provides an overview of AN/SPY-1 
specifications and variants installed on U.S. destroyers and cruisers as well as warships 
from five other countries. 


J J j j 


Attribute SPY-IB/D SPY-1D(V) SPY-IF SPY-1F(V) 


Primary 

Ship 

Classes 

Destroyers 

Cruisers 

Fngates 

Destroyers 

Cruisers 

Fngates 

Frigates 

Destroyers 

Amphibs 

earners 

Fngates 

Destroyers 

Amphibs 

Carriers 

Antenna 

Diameter 

12 ft. (3.7 m) 

12 ft. (3.7 m) 

8 ft. (2.4 m) 

8 ft. (2.4 m) 

No. of Elements 

4.350 

4.350 

1.856 

1.856 

SM-2 and 

ESSM Capable 

Yes 

Yes 

Yes 

Yes 

All-Weather 

Performance 

Yes 

Yes 

Yes 

Yes 


Command 

Guidanco 




o 

Terminal 

Guidance 


Illuminate 



Detect 



Figure 1.2: AN/SPY-1 variants (from Lockheed Martin, 2009). 


1.1.2 Tactical Data Links 

To develop and maintain a common tactical picture (CTP) among components of a battle 
group, networked platforms must be able to distribute information from each sensor in 
near real time. Tactical data links (TDLs) provide this capability. A TDL is an electronic 
connection between military platforms that allows for the transfer of digital information. 
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TDLs typically employ radio communications and transmit in the high frequeney (HF) or 
ultrahigh frequency (UHF) range. Standardized message formats allow for encrypted direct 
communication between eomputerized taetical data systems at high data rates (Hura et ah, 
2000, pp. 107-121). The two primary TDLs in use by North Atlantie Treaty Organization 
(NATO) forees are Link-11 and Link-16. 

Link-16 is the primary TDL for the distribution of situational awareness information. It 
uses UHF line-of-sight (LOS) transeeivers to establish a network capable of data rates up 
to 238 kilobits/see. Link-11 has a relatively low data rate of 1200 bits/see, but operates over 
both HF and UHF, allowing for eommunications beyond LOS (U.S. Army, 2000). Link-22 
is eurrently under development as a replacement for Link-11, providing greater data rates 
and a more stable conneetion over a seeure, beyond-LOS TDL. 

These TDLs allow radar track data such as position, heading, speed, and identification to 
be shared across platforms at rates up to 4,000 traeks per seeond (Jane’s, 2011). The traeks 
held by eaeh platform are compared to determine whieh traeks eorrespond to the same 
target. The radar tracking system automatieally eorrelates tracks that it determines are 
monitoring the same target. It then eombines the tracks and all available information into 
a single fused traek. This process has proven effeetive, but advances in proeessor speeds 
and TDL data rates provide the opportunity to aehieve a superior tactical picture through 
the Cooperative Engagement Capability (CEC). 


Cooperative Engagement Capability 

Rather than transmit processed track data, CEC allows distributed sensors to share unfil¬ 
tered measurements from each observation, such as range, bearing, and elevation. TDEs 
distribute these measurements throughout the battle group within mieroseeonds, allowing 
each unit to utilize the data as if it were provided by onboard sensors. With each platform 
performing the same traeking algorithm on all available data, an identical CTP is available 
to eaeh CEC platform (Johns Hopkins Applied Physies Eab, 1995). As Eigure 1.3 illus¬ 
trates, eombining radar measurement data to produce composite tracks results in a more 
aceurate contact picture than is possible from correlating multiple traeks aeross platforms. 
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Composite tracking 
and identification 



Figure 1.3: Composite tracking with CEC (from Johns Hopkins Applied Physics Lab, 1995, 
p. 379). 


1.2 Literature Review 

The problems associated with correlating multiple-target tracking across multiple sensors 
are well documented and analyzed. Bar-Shalom and Li (1995) provide an overview of 
many of the principles and techniques developed to address this topic. With the deployment 
of MLRs connected by high speed TDLs in the 1980s, the coordinated management of 
finite sensor resources to track an increasing number of targets became an active topic of 
research. 


1.2.1 Centralized Control of Radar Resource Management 

The majority of radar resource management literature focuses on the scheduling and tasking 
of a single sensor, such as in Bogler (1990) and Blac km an (1986). Many of the techniques 
developed are feasible in multiple-radar environments controlled by a central scheduler, 
though data sharing rates limit the feasibility of their implementation. Most early research 
on resource management for multiple radars utilizes a centralized control element to assign 
tasks to networked sensors based on task priorities and the sensors’ reported constraints 
and capabilities. Nash (1977) provides an early example of such approaches. 
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Given “a set of previously aequired and externally prioritized surveillanee targets,” Nash 
presents an algorithm for the optimal alloeation of targets to multiple sensors during the 
next time step. To do so, he utilizes the elassie transportation problem from linear program¬ 
ming (LP) which schedules the delivery of commodities from supply centers to customers 
at minimal cost to the system. In this case, the algorithm assigns sensors to targets in order 
to minimize the total cost incurred over the next time step. Costs are defined as a function 
of a target’s priority and the estimated track accuracy obtained from the sensor assignment. 
Track assignments are limited by the maximum track capacity of each sensor. Nash mod¬ 
ifies the basic transportation problem to allow more than one sensor to track a target. A 
benefit of using a version of the transportation problem is the availability of fast, heuristic 
approaches that allow for near real time calculations. Another benefit of this approach is 
the ability to modify sensor allocation constraints, such as by specifying target assignments 
or modeling reduced capacity. 

The algorithm provided by Nash is step-wise optimal when the assumptions of indepen¬ 
dence and linearity of track to sensor assignment costs hold. Unfortunately, this greedy 
approach does not imply optimality over the total time frame of observations. Since the 
formulation optimizes the gains obtained over the next time step vice all future time steps, 
a suboptimal allocation of tracking resources may occur. The author acknowledges this 
limitation, but suggests one can find a globally optimal solution by integrating the formu¬ 
lation within a dynamic programming process. 

Schmaedeke (1993) revisits this formulation, but calculates the predicted gains from 
sensor-to-target assignments with data available from the Kalman filter (KF). The KF is 
a recursive algorithm for estimating a target’s physical state based on a series of measure¬ 
ments over time. For each observation, the algorithm updates a multivariate distribution for 
each state variable being tracked. As a means of measuring the potential benefits of sensor 
assignments, the KF provides a measure of the information content in a target’s probability 
density function (PDF). Since it is possible to calculate the updated PDF a priori when us¬ 
ing a KF, the information content of a target’s PDF after the next time step can be calculated 
for each possible sensor assignment. Using the LP formulation developed by Nash (1977), 
one can determine the optimal sensor assignments. McIntyre and Hintz (1996) describe a 
search and track computer simulation that demonstrates the benefits of this approach over 
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random sensor-to-target assignments. 


By defining track accuracy estimates directly from the KF, Schmaedeke provides a rational 
and concrete means of determining the benefits provided by sensor to target assignments. 
However, suboptimal allocation of resources may still occur due to the formulation’s step¬ 
wise approach. Furthermore, defining the information content with the KF introduces ad¬ 
ditional restrictions. The functions developed do not provide for prioritizing targets or 
weighting a target’s state variables when deciding the value of information gained from 
sensor assignment. Additionally, for a traditional KF, the calculated information gain is 
unaffected by target maneuvers. To increase the value assigned to observations of ma¬ 
neuvering targets, Schmaedeke and Kastella (1998) recommend replacing the KF with in¬ 
teracting multiple model Kalman filters (IMMKF). By tracking and weighting several KF 
motion models continuously, the IMMKF provides a computationally intensive but more 
robust track estimate. During maneuvers, there is greater uncertainty as to which of the 
KF models is accurately tracking the target. This results in a reduction of the information 
content present in the IMMKF model. Thus, the possible information gain from a sensor 
assignment increases. This makes sensor assignment to maneuvering targets more valuable 
and thus more likely to occur. 

Rather than rely on the strict calculation of information gain prescribed by Schmaedeke 
and Kastella (1998), Blackman (1986, pp. 397-401) suggests employing utility theory 
to determine sensor allocation. To maximize overall utility, the algorithm calculates the 
marginal utility of assigning a sensor to a given track for each possible pairing. The util¬ 
ity of a target’s track depends on the value placed on different levels of track accuracy. 
The marginal utility is the difference in expected utility after the next time step with or 
without the sensor to target assignment. The system designer can control the shape of the 
utility function to ensure sensor assignments only add significantly to the marginal utility 
if they provide actionable level target accuracy. The marginal utility is also dependent on 
the weighted probabilities of a target’s type and threat assessment, both before and after 
sensor assignment. After calculating the marginal utility for each sensor-target pair, an op¬ 
timal assignment algorithm appoints sensors to targets. Unlike the approach presented by 
Nash, Blackman’s approach requires a one-to-one assignment of sensors to targets, which 
prohibits its direct application to real world tracking problems. 
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Veeramachaneni and Osadciw (2006) shift the focus from optimizing the next scheduled 
observation time to a more comprehensive approach for multiple sensor resource alloca¬ 
tion. Rather than optimizing a target’s track accuracy, they focus on the allocation of avail¬ 
able radar dwell time and energy between search and track functions for a sensor network. 
The objective function is the weighted sum of two competing objectives: the search cov¬ 
erage area and the probability of detection of tracked targets. Each time the target picture 
transforms, the controlling element must resolve the problem and reassign the radars. To 
approach an optimal solution in near real time, the authors suggest implementation of a 
particle swarm optimization algorithm (Eberhart & Kennedy, 1995). While this heuristic 
approach does not guarantee optimality, it does allow for continuously updated solutions to 
provide a best known allocation of time and energy in near real time situations. 

To overcome the limitations of step-wise optimal assignments, this thesis begins by assum¬ 
ing universal knowledge of targets. This includes knowledge of when the sensor network 
gains and loses a track and the target positions during this timeframe. This assumption 
provides the ability to determine an optimal observation schedule for each target that min¬ 
imizes the use of radar resources over the entire timeframe considered. While this is not 
a feasible approach for real radar scheduling, it does provide a “best case” bound on the 
resources required to maintain a given track picture. We relax this assumption later in 
the thesis when we investigate geographical partitioning methods. Unlike Veeramachaneni 
and Osadciw, this thesis assumes there is a specified amount of radar resources available 
for tracking. Eor simplification, we do not analyze the interrelationship of tracking capacity 
with other radar functions. 

1.2.2 Distributed Control of Radar Resource Management 

While centralized radar resource management techniques can provide optimal solutions, 
data sharing rates and computation speeds limit their feasibility in a distributed sensor 
network. Recent research attempts to overcome these technical limitations by develop¬ 
ing decentralized decision making algorithms. By developing distributed solutions that 
can approach optimality over time, the algorithms provide actionable decisions even when 
communications are limited. 

Dorman, Eeung, Nicholson, Siva, and Williams (2005) suggest two decentralized decision 
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making algorithms for sensor to target assignments and weapon allocations based on infor¬ 
mation theory. The first expands on the approach developed by Nash (1977). Instead of 
delivering all cost and threat data to a central planner to solve the LP problem and distribute 
target assignments, the network distributes the necessary data to all platforms which inde¬ 
pendently solve the LP to determine the targets for which they are responsible. While this 
approach can guarantee optimality within a tactical time frame for many networks, delays 
in data sharing can prevent the algorithm from reaching an optimal solution. 

To reduce the volume of data shared between platforms, Dormon et al. (2005) developed 
a second algorithm based on a decentralized auction. Rather than transfer all cost and 
threat data between each set of platforms, sensors only exchange one bid for each target. 
By limiting the information that must be shared, this algorithm reduces the possibility of 
exceeding a network’s data transfer capacity. The bids are used by each platform to run 
the decentralized auction algorithm independently and determine the optimal target assign¬ 
ments. The current framework for target data fusion across TDLs could allow these bids to 
be distributed by slightly reformatting the standardized target messages. This would limit 
the additional throughput capacity required for implementation and ensure the algorithm is 
scalable for situations with greater target density. 

Additional approaches to decentralized decision making rely on self-organization pro¬ 
cesses. Self-organization, or emergent behavior, is a process where meaningful global 
patterns emerge from the interactions of simple agents with each other and their local envi¬ 
ronment using only local information and without any centralized control (Camazine et ah, 
2001, pp. 7-8). For radar networks, the agents are sensor schedulers and the local environ¬ 
ment is the shared track picture. Stigmergy, a subset of self-organization in which agents 
only interact indirectly with one another through their impact on the environment, underlies 
the resource management approaches proposed in Weir and Sokol (2009) and Lambert and 
Sinno (2011). 

Weir and Sokol (2009) suggest local decision rules for radar management based on data 
from a shared track file. These rules result in global patterns consistent with optimal results 
obtained from centralized control. To determine which tracks should be observed over the 
upcoming time period, each platform independently develops scores for each target based 
on potential track accuracy improvements. The proposed scoring function is dependent on 
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the track state, track error covariance, range, and measurement uncertainties. Each sensor 
ranks the tracks to ensure the most effective observations happen first. Following these 
observations, the platforms update their shared track picture and score and rank the targets 
again. Appropriate local scoring criteria result in radar resource decisions that approach 
optimal solutions without central decision making or direct interaction between platforms. 

Beyond the potential to approach optimal resource allocation, Weir and Sokol (2010) also 
demonstrate that self-organizing radar managers react rationally to network degradations. 
When disruptions to network capabilities occur, the radars return to optimal decisions for 
independent operations. If one sensor loses communication with the network, the same 
local scoring criteria adapt to the new network configuration and radar resource decisions 
approach optimality over time. 

Rather than focus on a sensor’s potential impact on the future track picture, Lambert and 
Sinno (2011) develop local rules for resource allocation based on the information gain 
achieved for a given target during previous observations. Though radars do not interact di¬ 
rectly in stigmergic systems, they can sense the impact of others through their environment. 
By comparing the local information gain to the global information gain for an observed 
target, a radar can assess how important its contribution is to the global track picture. If 
the local information gain is significantly smaller than that observed in the global track file, 
Lambert and Sinno suggest measures to induce the radar to focus on other targets for which 
it may have a greater impact. 

These algorithms developed for distributed control of resources still require a significant 
amount of data transfer and computational effort each time sensor-to-target assignments 
are optimized. A portion of this thesis investigates the effect of assigning sensors to targets 
based on geographic location. The minimal computational requirements of this method 
are conducive to real time implementation. However, reliance on geographic partitioning 
alone oversimplifies the problem, leading to periods where sensors lack the resources to 
adequately monitor their assigned targets. While not a replacement for the approaches 
mentioned in this section, the assignment methods developed in this thesis may be able to 
reduce the real time computational demands of these algorithms by reducing the problem 
size. 
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1.3 Problem Statement and Assumptions 


This thesis investigates how to manage a distributed network of radar assets when operating 
in a dense target environment. Several factors complicate the optimal use of available radar 
resources. The multifunctional nature of modem radars places demands on resources often 
not known to the network. Furthermore, multiple objectives exist when scheduling radar 
resources, namely minimizing emitted power and dwell time. Optimization of one of these 
does not necessarily result in optimal or near optimal results for the other. To continue 
exploring this problem, we focus solely on the tracking function of MFRs and make some 
simplifying assumptions. 

To remove directional and environmental effects, we assume variations in sensor effective¬ 
ness only depend on range and emitted power. That is, we do not consider the complex 
effects of weather, land objects, and the directional limitations of MFR arrays on radar 
tracking. We also assume that all targets have identical tracking requirements; that is, 
we do not attempt to utilize knowledge of targets to prioritize them by threat levels. The 
sensor-to-target assignment models in this thesis do not rely on these assumptions, but the 
scenarios considered in our analysis comply with them. These assumptions are necessary 
for the partitioning models presented. 

Common to all models, we assume that any radar return above a signal to noise ratio (SNR) 
threshold of 13 dB provides the same information gain to the sensor network. This means 
that there are no false positive or false negative returns, that the radar assigns each return 
correctly to its track, and that the resultant track quality after any observation is equivalent. 
Finally, we assume that track data is distributed to all sensors by the TDL without time 
delay and unrestricted by throughput capacity. 

This thesis is organized as follows. Chapter 2 introduces an optimal radar to target assign¬ 
ment model that provides a best-case radar schedule for networked MFRs. We extend this 
model to allow for analysis of the potential difficulties introduced to the radar scheduling 
process by the coordinated maneuvers of targets. Chapter 3 introduces an optimal partition¬ 
ing model for assigning tracking requirements to sensors based on the geographic location 
of targets. We also present two variations to this model. The first permits regions of space 
to remain unassigned, allowing any sensor to perform tracking in these regions. The sec- 
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ond accounts for the potential difficulties of handing off targets between sensors. Chapter 
4 presents the analysis performed using each of these models, and Chapter 5 provides a 
conclusion and recommendations for further research. 
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CHAPTER 2: 

OPTIMAL RADAR ASSIGNMENT IN A KNOWN 
TARGET ENVIRONMENT 


In this chapter, we investigate the performance bounds of networked sensors. Specifically, 
we develop a model to determine the best use of radar resources when multiple MFRs 
coordinate through an infinite-capacity network with zero time delay. We consider an envi¬ 
ronment in which the sensor scheduler has universal knowledge of each target’s locations 
and maneuvers over the period of time tracking occurs. Eliminating the uncertainty in fu¬ 
ture events allows the scheduler to optimally allocate resources over the entire period of 
interest. While this algorithm is not deployable as a real-world radar scheduling technique, 
it does provide a best-case bound to which one can compare other algorithms. 

2.1 Multifunction Radar Capabilities 

For a radar to recognize a return from a target, the signal to noise ratio (SNR) of the tar¬ 
get echo must exceed an established threshold. A number of factors effect how much of 
a radar’s emitted energy is sensed by the radar receiver. External variables, such as the 
target’s range and environmental noise, as well as internal factors, such as antenna gain 
and electrical losses, effect how strong the radar signal is compared to all sources of noise. 
Equation (2.1) is a reproduction of the radar equation provided in Sabatini and Tarantino 
(1994, p. 27). Each factor presented, whether controlled by the radar or not, can impact the 
SNR. 


SNR 


\trans ^trails ^rec ^2 ptrans prec2 




( 2 . 1 ) 


where 

ptrans 

pprans 

Qvec 


Transmitted power [MW] 

Transmitter antenna gain in the direction of the target [dimensionless] 
Receiver antenna gain in the direction of the target [dimensionless] 
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A 

ptrans 

prec 

a" 

R 

kg 

To 

pn 

ptot 

T 


Transmit wavelength [em] 

Transmitter propagation factor [dimensionless] 
Receiver propagation factor [dimensionless] 
Target radar cross section (RCS) [m^] 

Range to the target [m] 

Boltzmann constant [1.38vlO~^^J/K] 

Standard Temperature [290 K] 

Receiver noise figure [dimensionless] 

Total losses for the radar system [dimensionless] 
Radar pulse width [/isec] 


Table 2.1 provides values for the variables in Equation (2.1) derived from the benchmark 
problem for phased array radars developed in Blair and Watson (1996). The radar cross 
section (RCS) of all targets is established as one square meter for simplification, but the 
following formulations are not limited to this value or to a constant RCS throughout the 
trajectory of the targets. Blair and Watson establish eight discrete waveforms available for 
use with varying radar pulse widths, tk- 


ptrans 

1 MW 


No. 

T (/isec) 

Qtrans 

4752 


1 

1.35 

Qrec 

4752 


2 

2.25 

A 

7.5 cm 


3 

4.05 

ptrans 

1 


4 

5.85 

prec 

1 


5 

11.70 

o° 

1 m 


6 

23.40 

pn 

2 


7 

46.8 

ptot 

144.5 


8 

93.6 


Table 2.1: Radar parameters and pulse width options (after Blair and Watson 1996). 


The radar pulse width must be long enough to ensure adequate energy is emitted to obtain 
an SNR above the sensor’s threshold. We assume that an SNR above 13 dB provides perfect 
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detection without false alarms. Furthermore, we assume that any SNR above the threshold 
results in the same track accuracy and requires the same update rate. Thus, the radar utili z es 
the waveform that meets the SNR threshold with the least emitted energy. 

2.2 Sensor-to-Target Assignment Model 

We now describe a model that optimizes the schedule of observations required from net¬ 
worked sensors for a given set of tracks. For this model, we make the additional assumption 
that we have prior knowledge of each track’s duration and location. This allows the model 
to determine the optimal schedule over the entire time horizon. To limit computational 
effort, we establish a dynamic set of waveforms from each sensor that are available for 
observing each target at each time step. Since we assume each observation provides the 
same information gain, the set only contains the waveform from each sensor that meets the 
SNR threshold with the least emitted energy. 

Decision variables describe the optimal scheduling of track observations by sensors 
to minimize the amount of power and time invested in tracking and limit the impact of 
delays due to radar resource limitations. We assume that longer time delays result in a 
continued reduction in track quality and an increasing probability of track loss. The model 
utilizes multiple linear constraints to capture this nonlinear impact on the quality of a track¬ 
ing schedule. This ensures the model limits both the total time of all delays and the length 
of those delays. 

2.2.1 Formulation 

Index Use rCardinality] 


t G T 

time steps [~100] 

i G I 

sensors ["3] 

j^J 

targets [~30] 

keK 

waveforms [~8] 

1 G L 

time steps during the required observation interval L [~10] 

neN 

segments of a piece-wise linear delay penalty function[~2] 

{j,t) eM cJxT 

time steps t in which target j is tracked ["3000] 
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Given Data [Units] 


Ptransmit 

interleave_limit 

nin 


pulse width of waveform k [/isec] 
transmitter power [MW] 

number of targets a sensor can observe per time step [targets] 
slope for penalty constraint n [J] 
intercept for penalty constraint n [J] 


Calculated Data ['Cardinality] 

(z, j^t^k) e P C I X M X K set of waveforms k from sensor i that are capable of ob¬ 
taining an SNR above the threshold SNRmin for target j at 
time t with the least emitted energy [~9000] 


P= { {iJ,t,k)\{j,t) GMand 


k = argmin 

k 




ptransQtransQrec ^2ptransprec2 


i,J,l 


CT, 


i,J,l i,J,t i,J,t 


(4,,)Xuk«ToF,"L'% 


Tk > SNR„ 


Decision Variables [Units] 

binary decision to observe target j with sensor i at time t with 
waveform k [binary] 

NO_LOOKj t binary decision to miss required observation of target j between 

times t -I- 1 and t + \L\ [binary] 

DELAYj t number of time steps by which an observation of target j has been 

delayed as of time t -I- |L| -|- 1 according to the observation fre¬ 
quency requirement [time steps] 

DELAY_PENALTYj t penalty for number of time steps that an observation of target j at 

time t -|- |L| -|-1 is delayed beyond the required observation interval 
\L\ [J] 
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Formulation 


min £ [ptransmitiTk-j)Xij,,,k)+ £ DELAY_PENALTYj,t (2.2) 

NO_LdOK, {j,t)eM 

DELAY, 

DELAY _PENALTY 


S.t. 


£ < interleave_limit 

j,k:{i,j,t,k)eP 

£ Xij^t+i,k+NO_LOOKj,t>( 

l,i,k:{i,j,t,k)eP 

DELAYj^t-DELAYjj I > 1 - (1 - NO_LOOKj^t)t 
DELAY_PENALTYj^t > mn(DELAYj^t - DELAYj^t+i) 
N0_L00Kj^t, delaydelay_PENALTY j^t > 0 

e {0.1} 



'ii,t 

(2.3) 

V(7T) 

eM 

(2.4) 

v(jt) 

eM 

(2.5) 


eM 

(2.6) 

v(jt) 

eM 

(2.7) 


\eP 

(2.8) 


2.2.2 Discussion 

The objective function (2.2) expresses the total cost of a solution as the sum of the energy 
transmitted and the penalties for delaying observations of a target. To eliminate multi¬ 
ple optimal solutions and reduce computational time, the objective includes the factor y to 
force sensor observations to occur as early as possible without requiring additional observa¬ 
tions or more costly waveforms. DELAY_PENALTYj t models a piecewise linear, convex 
penalty function that is large enough to ensure an observation is performed, if possible, 
rather than accept the penalty. 

Each constraint (2.3) prevents the scheduling of observations by a single sensor at a given 
time in excess of the interleaving limit. To ensure each target is monitored at the required 
sampling frequency, each constraint (2.4) requires at least one observation to occur within 
\L\ time steps of t or records the missed observation. 

Constraints (2.5) and (2.6) combine to calculate the penalty for delaying observations. Con¬ 
straints (2.5) determine the amount of time by which observations of target j are delayed 
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at time t+ \L\. If NO_LOOKj t is zero, then DELAYj t is zero. Otherwise, DELAYj t is 
the number of consecutive time steps NO_LOOKj t has been equal to one up to time t. 
Constraints (2.6) utilize DELAY jt to calculate the penalty for delaying observations. The 
constraints only assess a penalty for the time that a delayed observation actually occurs. For 
those time steps with a delayed observation, the difference between consecutive DELAY jt 
is equivalent to the number of time steps that the schedule delays the observation. As 
Figure 2.1 illustrates, multiple linear constraints combine to obtain a non-linear penalty 
response. This causes one long delay to be more expensive than multiple shorter delays 
of the same total length. Table 2.2 provides an example of the calculations performed by 
constraints (2.5) and (2.6). Constraints (2.7) and (2.8) define the domains of the decision 
variables. 


DELAY PENALTY Constraints 



Figure 2.1: Delay penalty constraints. 


t 

... 14 

15 

16 

17 

18 

19 ■■■ 

NO_LOOKjj 

... 0 

0 

1 

1 

1 

0 ■■■ 

DELAYjj 

... 0 

0 

1 

2 

3 

0 ■■■ 

DELAY _PENALTYj^t 

... 0 

0 

0 

0 

5000 

0 ■■■ 


Table 2.2: Notional delay penalty calculation. 
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Figure 2.2 provides a geographical representation of results from the sensor-to-target as¬ 
signment model for a network of three sensors maintaining 20 tracks. The sensors, rep¬ 
resented by circles, are stationary throughout the problem and the targets, represented by 
black arrows, move along their predetermined tracks at constant speed and direction. The 
hash marks are observations along the targets’ paths. The size and color of the hash marks 
correspond to the power utilized and the sensor performing the observation, respectively. 



Figure 2.2: Geographic display of results for scenario involving three sensors and 20 tar¬ 
gets. 


Figure 2.3 (a) provides a timeline of the results presented in Figure 2.2. For this example, 
the interleave_limit is set to one, constraining each sensor to observe at most one target 
per time period. The observation frequency \L\ is 10 time steps. Figure 2.3 (b) displays 
the power requirements for each sensor to execute the optimal observation schedule as a 
running average over 10 time steps. 
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Figure 2.3: Timeline and power demands of optimal results for scenario involving three 
sensors and 20 targets. 


2.3 Maneuverable Targets 

To investigate the potential effects of airborne adversaries capable of coordinating their 
operations, we extend the sensor-to-target assignment model to allow for target maneu¬ 
vers. Following a target maneuver, the sensor network must perform additional observa¬ 
tions to account for the reduction in track accuracy. To model this, we declare the subset 
(7: t) & QcMio define the time t that target j maneuvers. We then increase the number of 
observations required during the period of length \L\ immediately following the maneuver. 
For each of the observations following a maneuver that the sensor network is unable to 
schedule, the formulation assesses a penalty. 
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2.3.1 Formulation 

In addition to the indices, data, and variables described in Section 2.2, we introduce the 
following nomenclature. 

Index Use [~Cardinality] 

{i,t) ^ QC M time step t at which target j maneuvers [“30] 

Given Data [Units] 

turn_looks number of observations required during normal observation fre¬ 
quency \L\ following a target maneuver [observations] 
turn_penalty penalty for number of observations missed following a target ma¬ 
neuver [J/observation] 

Decision Variables [Units] 

MISSED_LOOKSj number of required observations missed follow¬ 
ing the maneuver of target j [observations] 


Formulation 


Y Arr^Yhn)<r ^ iptransmiti'^k + 52 DELAY_PENALTYj t 

DELAY. ’ 

DELAY _PENALLY, 

MISSED_LOOKSj 


+ turn_penalty^MISSED_LOOKS j 
j 


(2.9) 
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s.t. 


Xi,j,t+l,k + MISSED_LOOKSj >tum_looks 

{l,i,k,t):{j,t)eQ and {i,j,t,k)eP 


I 


^i,j,t,k 


{j,k):{i,j,t,k)eP 


< interleave 


^ Xij^t+i^k + NO_LOOKj,t>\ 

{l,i,k):{iJ,t,k)eP 


DELAYj^t- delay > 1 - (1 -NO_LOOKj^t)t 
DELAY_PENALTYj^t > mn{DELAYj^t -DELAYj^t+i) +K 


NO_LOOKj^t:DELAYj^t:DELAY_PENALTYj^t > 0 


MISSED_LOOKSj > 0 


Xi,j,t,ke{OA} 


Vj 

( 2 . 10 ) 

\/iA 

( 2 . 11 ) 

y{j,t)eM 

( 2 . 12 ) 

y{j,t)eM 

(2.13) 

\/n,{jA) eM 

(2.14) 

y{j,t)eM 

(2.15) 

y{j,t)eM 

(2.16) 

V/, (;V) eM,k 

(2.17) 


2.3.2 Discussion 


The objective function (2.9) expresses the total cost of a solution as the sum of the en¬ 
ergy transmitted and the penalties for delaying or missing observations on a target. It is 
similar to the objective function (2.2) from Section 2.1 with one additional component. 

turn_penaltyY^MISSED_LOOKSj penalizes a solution for failing to perform the addi- 
j 

tional observations required following a target’s maneuver. As with DELAY_PENALLYj t, 
the turn_penalty is large enough to ensure observations are performed, if feasible, rather 
than accept the penalty. 
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Constraints (2.10) reflect the required turn_looks during the period \L\ following a ma¬ 
neuver by target j at time t. Each additional observation that the set of sensors fails to 
schedule during this timeframe is tallied as MISSED_LOOKSj. Constraints (2.11) prevent 
the scheduling of observations by a single sensor at a given time in excess of the interleav¬ 
ing limit. To ensure each target is monitored at the required non-maneuvering sampling 
frequency, each constraint (2.12) requires at least one observation to occur within \L\ time 
steps of t or record the missed observation. 

Constraints (2.13) and (2.14) combine to calculate the penalty for delaying observations 
and are identical to constraints (2.5) and (2.6) discussed in Section 2.2.2. Constraints 
(2.15),(2.16) and (2.17) define the domains of the decision variables. 

The computational requirements of this problem preclude an exhaustive analysis of net¬ 
works of significant size or complexity. In an attempt to reduce the computational effort, 
we developed an attacker-defender model that allowed for one maneuver per target over the 
time horizon. The formulation failed to provide computational benefits because the solver 
had to execute the subproblem for each possible maneuver schedule before obtaining op¬ 
timal results. To remove the unnecessary calculations of maneuver schemes, we instead 
iterated through every feasible combination of maneuvers for the targets to determine the 
“worst case” set of maneuvers. 

For simplification, we maintain the assumption that each target may perform at most one 
maneuver during the planning horizon. Furthermore, for each time step, we assume there is 
only one maneuver available to each target being tracked. For example, at any time step, a 
target only has the option to turn 90 to the right. Figure 2.4 displays the schedule obtained 
following the “worst case” maneuvers for a scenario with two sensors and six targets arriv¬ 
ing every two time steps. The time horizon for this scenario is 20 time steps. The sensors 
must observe each target at least once every three time steps. Following a target maneu¬ 
ver, three observations must occur within the next three time steps or a penalty is assessed. 
We allow for target maneuvers any time from one time step after target tracking begins 
to three time steps before the target track terminates. By preventing later maneuvers, we 
ensure that the sensors have the opportunity to perform the required observations following 
a target maneuver. In the “worst case” maneuver scheme, the targets synchronize their ma¬ 
neuvers into two blocks to concentrate the tracking demands on the sensor network. The 
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decision of the targets to maneuver in different time blocks is a result of the varying track 
arrival and departure times. In this scenario, the coordinated decisions of targets result 
in insufficient sensor capacity to perform all required observations following maneuvers. 
Section 4.1 discusses additional results from this formulation. 



Time Step 


Figure 2.4: Optimal observation schedule with “worst case” maneuvers for scenario in¬ 
volving two sensors and six targets. 
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CHAPTER 3: 

TARGET ALLOCATION BASED ON GEOGRAPHIC 

LOCATION 


The formulations of Chapter 2 rely on the assumption that the sensors have perfeet knowl¬ 
edge of the target locations and track durations. This is useful for identifying trends in 
the resultant schedules but prohibits the application of the model in real world situations. 
Interestingly, results obtained from these formulations suggest significant regions of the 
target operating space are often dominated by a single sensor. Figure 3.1 displays the re¬ 
sults of 120 simulations geographically. Each simulation distributes 25 targets randomly 
throughout the operating space while the locations and capabilities of the three sensors re¬ 
main constant. Observations indicate that a single sensor begins to dominate regions as the 
distance to the sensor network increases. This suggests that a target’s physical position may 
provide enough information for a network to allocate some portion of targets to sensors. 




Figure 3.1: Observations from optimal schedules for 120 simulations involving three sen¬ 
sors and 25 targets. 
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To investigate the possibility of assigning radars to targets based on geographic location, we 
develop a model to partition the physical decision space so as to balance the resultant radar 
resource demands between sensors. This model does not rely on the assumption that sen¬ 
sors have perfect knowledge of the targets. Instead, it relies on an estimate of the expected 
target density within the operating space. To account for the increased uncertainty intro¬ 
duced by eliminating this assumption, we suggest a relaxation to this partitioning model 
that limits the amount of space assigned to a specific sensor and allows any sensor to track 
targets in the unassigned regions. This reduces the probability of exceeding a sensor’s 
tracking capacity by limiting the sensor-to-target assignments prescribed based on geo¬ 
graphic location. To investigate the impact that transferring target tracking requirements 
between sensors has on geographical partitioning, we develop another model to account 
for the increased workload inherent in target handoffs. 


3.1 Geographic Assignment 

To analyze the potential benefits of assigning sensors to targets geographically, we extend 
an optimal partitioning model developed by Carlsson (2012) to three dimensions. In addi¬ 
tion, we adjust the formulation to account for radar power requirements, individual sensor 
tracking capacity, and varying target densities throughout the physical decision space. 

3.1.1 Formulation 

To partition space between sensors, we consider three factors. First, to observe a target 
at position x, a sensor must emit energy proportional to the fourth power of the range to 
the target, as required by Equation (2.1). Second, the assigned regions of space should be 
proportional to the capacity that each sensor can provide to tracking functions. Radars in 
a degraded status or involved heavily in other functions, such as searching or illuminat¬ 
ing, should handle less of the tracking demands for the sensor network. Third, points in 
space should be weighted by the expectation of observing a target at that position during 
a scenario. Regions of high target density will result in proportionally greater demands on 
radar resources than areas of low target activity. The following mathematical formulation 
accounts for these factors. 
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Index Use TCardinality] 


jc G Z) positions in physical space [infinite] 
i G I sensors [~3] 

Given Data [Units] 

at relative eapaeity of sensor i to eonduet traeking 

[energy/target/dist^] 

p{x) expeeted density of targets evaluated at position x [tar¬ 
gets/time/volume] 

Pi position of sensor i [unitless] 

Decision Variables [Units] 

Ri region of spaee assigned to sensor i [unitless] 


Formulation 


min max 

^ 1 r • • ^ 



aip{x)\\x-pi\\‘^dV 


(3.1) 


s.t. 

n 

= ^ (3-2) 

/=1 


3.1.2 Discussion 

The objeetive funetion (3.1) minimizes the maximum demand on any sensor by ehoosing 
whieh regions of spaee to assign to eaeh sensor. The integrand aeeounts for the sensor 
and target attributes previously mentioned as eritieal to aequiring an optimal partition. The 
density funetion p{x) can aeeount for different target types by utilizing Equation (3.3). 
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(3.3) 


p(^) = E Pjpjix) 

j=i 

Parameter j8; is a weighting for target type j to account for unique factors, such as RCS, 
threat level, and update requirements. The constraint (3.2) ensures all of the physical space 
is assigned to at least one of the sensors’ regions. 

3.1.3 Dual Formulation 

The formulation of Section 3.1.1 permits reformulation as an infinite-dimensional mixed 
integer linear programming (MILP) problem. The following formulation is the dual of 
the LP problem obtained by relaxing the integer constraint of the MILP problem. This 
formulation provides a means of computing the solution to the formulation presented in 
Section 3.1.1 by reducing the number of decision variables to a cardinality proportional to 
the number of sensors. Carlsson (2012) provides the proof of duality, which is reproduced 
in the appendix, and a computational algorithm. Indices and given data remain consistent 
with Section 3.1.1. 


Decision Variables [Units] 

Xi dual variable associated with sensor i [unitless] 

o{x) dual variable for tracking sensor assignment to position x [m^] 


Formulation 


max 
A,a 



D 


(3.4) 


s.t. 
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o{x) < ?iiaip{x)\\x-pi\\^ 

\/x, i 

(3.5) 

VI 

<-< 


(3.6) 

A/>0 

Vz 

(3.7) 


3.1.4 Dual Discussion 

The dual objective function (3.4) maximizes the sum of sensor assignment costs across all 
space D. Constraints (3.5) bound o at position v: by the assignment costs for each sensor. 
Constraints (3.6) ensure the sum of all A, do not exceed one. Constraints (3.7) define the 
domains of the decision variables A/. 

Assigning regions with Equation (3.8) results in optimal partitioning. The appendix pro¬ 
vides an adaptation of the theorem and proof from Carlsson (2012) for our formulation. 
Each location v: in D is assigned to sensor i only if the cost A, a,j|a: —p,j|^ is less than the 
cost of assignment to any other sensor. 


Ri = ^xGD I A;a,j|jc —p,j|^ < Ajay||.r — Vj7^z| (3.8) 

This formulation provides several benefits for a distributed radar network. Eirst, it permits 
inclusion of the critical factors established in Section 3.1.1. Second, it is scalable for the 
inclusion and removal of sensors. The dynamic nature of radar networks in an adverse 
environment requires the partitioning process to allow for changes in network size and 
geometry. This formulation allows for the recalculation of partitions at speeds consistent 
with the transitions of sensor network structures. Einally, the partitioning is adaptable to 
fluctuations in individual radar tracking capacity. If new demands or limitations on sensor i 
reduce the amount of resources available for tracking, the region assigned to that radar can 
be reduced proportionately by varying the value of a,-. A sensor’s declaration of reduced 
capacity through the TDL, or a loss of communication with the network, will allow the 
other sensors to calculate their newly assigned regions. 
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Figure 3.2 displays an expected target density function p(x) for a three dimensional sce¬ 
nario with a single HVA and three directions of approach for targets. Three networked 
sensors are positioned between the HVA and the threat axes. Applying the geographic 
partitioning model, each sensor is assigned the regions displayed in Figure 3.3. 




Figure 3.2: Three dimensional target density for one HVA with three directions of ap¬ 
proach. 




Figure 3.3: Optimal partition of space for target density in Figure 3.2. 
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3.2 Relaxation of Geographic Assignment 

By partitioning physical space for each sensor, the sensor-to-target assignment problem 
reduces to that of determining which region the target falls in. This provides for a signif¬ 
icant reduction of real time computational efforts for radar scheduling, but it also results 
in time periods that overwhelm individual sensor tracking capacities while under-utilizing 
other sensors. This results from intervals of time during which targets concentrate within a 
region at levels in excess of the time-averaged target density. 

To account for the time variance in target densities, we reduce the amount of space D as¬ 
signed to a specific sensor by including parameter yin Equation (3.8), resulting in Equation 
(3.9). Eocation v: in D is assigned to a partitioned region only if there exists a sensor i whose 
cost is less than the cost of assignment to all other sensors by at least y. This 

results in an unassigned region in which targets may be tracked by any sensor. 

Ri=^^xeD I ^iai\\x-piW"^ < XjttjWx-pj\\^-y VjV*} (3-9) 

Eigure 3.4 displays the results of a uniform average target density with varying values of 
y. As the magnitude of y increases, the inequality becomes more difficult to satisfy and 
Equation (3.9) assigns less of the original space D to a specific sensor. As the range from 
the sensor network increases, the inequality becomes easier to satisfy as the magnitude of 
Xiai\\x — pi\\‘^ becomes much greater than y. We present further results and analysis from 
this formulation in Section 4.2 
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Figure 3.4: Effects of yon the geographic partitioning results for three sensors in a uniform 
average target density. 


3.3 Geographic Assignment with Handoff Penalties 

In real world applications, transferring tracking responsibilities between sensors can con¬ 
sume additional radar resources. The uncertainty of a target’s relative bearing and range, 
the dependence of RCS on the target and sensor geometry, and the potential for environ¬ 
mental interference can require additional effort by the newly assigned radar to establish 
the required track quality. To include these additional radar resource costs in a partitioning 
scheme, we introduce a model that incorporates handoff penalties. 


34 













3.3.1 Formulation 

To account for handoffs, we discretize the formulation introdueed in Seetion 3.1.1 and 
add a penalty for expeeted target traeks entering a region. To assess the penalty, we define 
speeifie direetional target paths through the diseretized spaee. The model penalizes a sensor 
if one of these paths enters its assigned region. 

Without additional eonstraints, a target’s assigned region may inelude the speeified target 
paths to prevent a handoff penalty while ignoring the area between the path and the sensor. 
To prevent this, we enforee star eonvexity of the regions. For a diserete eell to be part of 
a sensor’s assigned region, all of the eells between it and the sensor must also be in that 
region. The resultant mathematieal formulation follows. 

Index Use [~Cardinality] 

i G I sensors [~3] 
j eJ targets [~30] 

k e K geographie position (alias k') [~2500] 

Sets 

(k, kf, j) e P k' is the diserete position that follows k along the path of target j 
(k, k'^ i) G N k' is the neighboring diserete position elosest to sensor i for position k 

Given Data [Units] 

Pi^k energy eost to monitor diserete position k with sensor i [J] 
pk target density at geographie position k [targets/s] 
a penalty for permitting handoffs [kW] 

/i weighting faetor to prioritize objeetives [unitless] 

Decision Variables [Units] 

Z energy demands on maximally tasked sensor [J] 

Xi^i^ binary deeision to assign position k to sensor i [binary] 

U, yt binary deeision to gain target j with sensor i at position k [binary] 
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Formulation 


M I y^. PkPi.k^i.k + CC ^ ^i,j,k I (3.10) 

’ ’ \i,k ij,k J 


^ PkPi,k^i,k T ^ 

'Y^yij,k 

Vz 

(3.11) 

k 

YiJ,k > Xi,k ~ ^i,k' 

j,k 

Vz, {k,k'j)eP 

(3.12) 

O 

VI 

1 


Vz, [kjk'j) G N 

(3.13) 

s 1 

i 


Vzt 

(3.14) 

yi,j,k,Xi,ke{0A} 


yij,k 

(3.15) 


3.3.2 Discussion 

The objective function (3.10) minimizes the energy demands on the maximally tasked sen¬ 
sor and, subordinately, minimizes the overall energy demand. Constraints (3.11) ensure 
that Z tracks the maximum energy demanded of any sensor by bounding it by the energy 
demands of tracking and handoffs for each sensor. To enforce a handoff penalty, constraints 
(3.12) require Yij^^ to be one if sensor i is assigned target j at position X, ^ but was not as¬ 
signed the previous location along the target’s path, Constraints (3.13) enforce star 
convexity of the assigned regions by ensuring that a discrete position k is not assigned to 
sensor i unless the path between i and k is also assigned to i. Constraints (3.14) ensure that 
each position is assigned to at least one sensor. Constraints (3.15) define the domains 
of the decision variables. 

To demonstrate the results of this formulation, we project the target density function in 
Figure 3.2 onto two dimensions and discretize the results. Figure 3.5 displays the new 
expected target density within a 50 x 50 hexagonal grid. The three black lines represent the 
target paths for which the model assesses penalties if a handoff between sensors occurs. 

Figure 3.6 presents the optimal results from this formulation when there is no handoff 
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penalty and when the penalty for handoffs is large. The results with no handoff penalty 
are similar to those we obtain with the continuous formulation of Section 3.1. When the 
penalty becomes adequately large, the cost of handoffs exceeds the costs of reassigning 
space to maintain track paths within a single sensor’s region. 



Figure 3.5: Two dimensional target density for one HVA with three directions of approach. 


No Handoff Penalty Large Handoff Penalty 



Figure 3.6: Results of geographic assignment with handoff penalties for the target density 
in Figure 3.5. 
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CHAPTER 4: 

RESULTS AND ANALYSIS 


This chapter presents the results and insights from formulations presented in Chapters 2 
and 3. We performed all computations on a Dell Precision T1650 PC with a 3.40 GHz 
Intel Core 17-3770 CPU and 16 GB of RAM. To solve the geographic assignment model of 
Section 3.1, we utilized MATLAB rev. R2012b. We implemented the other optimization 
models using the General Algebraic Modeling System (GAMS) rev. 236 and solved them 
with CPLEX 12.2.0.2. The details of specific problem instances, including approximate 
model sizes and computation times, are included in the following sections. 

4.1 Maneuverable Targets 

We investigate the impact of coordinated target maneuvers by using the model described in 
Section 2.3. For each scenario, the model contained approximately 250 decision variables, 
of which approximately 170 were discrete, and 200 constraints. Solution times were less 
than one second. To determine the “worst case” maneuver scheme, we considered 4096 
possible maneuver combinations. The total model run time averaged two hours. 

Figure 4.1 displays the schedules obtained for four scenarios. In each scenario, two sensors 
are tasked with tracking six targets. The sensors are capable of observing at most one 
target per time step. Upon commencement of tracking, the sensor network must observe 
each target at least once every three time steps. Thus, when tracking all six targets, the two 
sensors are operating at their maximum capacity. The sensors track each target for eight 
time steps. Following a target maneuver, three observations must occur within the next 
three time steps or a penalty is assessed. Maneuvers are not allowed after the fifth time step 
during which the target is tracked. This ensures the sensors have three time steps in which to 
perform the required observations. Each scenario varies only in the times at which tracking 
of the targets commence. In the first scenario, the sensors commence tracking all targets 
during the first time step. Target tracks commence at one, two and three second intervals 
in the second, third, and fourth scenario, respectively. For each scenario, we consider all 
feasible combinations of target maneuvers to determine which has the greatest impact on 
the radar tracking schedule. 


39 




Arrival Interval of Targets: Zero Time Step 



Sensor 1 
Sensor 2 
Turn 

Completed Turn Observations 
Missed Turn Observations 


Missed Turn Observations: 12 


Time step 

Arrival Interval of Targets: One Time Step 




Arrival Interval of Targets: Three Time Step 


« 3 
Ut 

1 



Missed Turn Observations: 0 


Time Step 


Figure 4.1: Optimal observation sehedules with “worst ease” maneuvers for four seenarios 
involving two sensors and six targets. 
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As demonstrated in the first seenario, the targets attempt to synehronize their maneuvers in 
order to maximize the demands for sensor assets in the following time steps. This results in 
a failure to perform suffieient observations on any target following the maneuver. As the in¬ 
terarrival time of target traeks inereases, the targets beeome less effeetive at overwhelming 
sensor resourees due to the inability to synehronize all of their maneuvers. By the fourth 
seenario, in whieh target traeks are offset at three seeond intervals, the targets are no longer 
able to eoordinate their maneuvers so as to overwhelm the sensors. 

These results suggest that adversarial targets eould overwhelm a sensor network operating 
at or near its traeking eapaeity by effeetively eoordinating their operations. The impaet on 
track quality of delayed observations, especially during this period of varying target speeds 
and directions, could result in the loss of a large number of hostile tracks. However, the 
effects diminish as targets lose the capability to effectively synchronize maneuvers. 


4.2 Geographic Assignment 

We now use the models from Sections 3.1 and 3.2 to analyze the performance of sensor- 
to-target assignments based on geographic partitioning. Each instance contained between 
200,000 and 450,000 decision variables, of which 5,000 to 15,000 were discrete. The num¬ 
ber of constraints ranged from 6,500 to 13,000. Typical solution times ranged from 10 
seconds to 30 minutes per simulation. Solution times increased significantly when obser¬ 
vation delays occurred. 

Figure 4.2 and Figure 4.3 present the results of six scenarios involving three stationary 
sensors. Each scenario consists of a different number of targets, from 25 to 33, that have 
randomly assigned arrival and departure times within the 105 time steps considered. A 
sensor can perform at most one observation per time step, and the sensor network must 
observe each target at least once every 10 time steps or accept a delay penalty. Thus, a delay 
occurs any time the sensor network is monitoring more than 30 targets simultaneously. 
We conduct 120 simulations of target tracks for each scenario. For each simulation, we 
determine the “best case” schedule of observations by implementing the sensor-to-target 
assignment model of Section 2.2. 
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Figure 4.2: Average delays, observations, and transmitted energy for scenarios involving 
various target densities. 
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Figure 4.3: Delays, observations, and transmitted energy for each simulation of scenarios 
involving various target densities. 
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To investigate the effeetiveness of the geographieal partitioning approaeh developed in See- 
tion 3.2, we choose various values for 7 and then solve for the “best case” schedule. The 
model restricts sensor assignments in the partitioned regions by only assigning waveforms 
from the designated sensor to the set of available observations (/, 7 , t, k) G P. For each value 
of 7 , we obtain the “best case” schedule for 120 simulations of each of the six scenarios. 

Figure 4.2 displays the average number of delays, average number of observations required 
above optimal, and the average amount of energy required above optimal for each value of 
7 and number of targets. Figure 4.3 displays the number of delays, number of observations 
above optimal, and amount of energy above optimal required for each simulation. For all 
scenarios, the impact on the “best case” schedule appears to be relatively stable until the 
partitioning method assigns approximately 40 percent of the target tracks to sensors. After 
this point, the number of delays increases significantly as targets overwhelm the sensors’ 
tracking capacities when they concentrate within a single sensor’s region of responsibility. 

We perform similar analyses to investigate the impact of different sensor network sizes and 
geometries. The sensor capabilities remain the same as in the previous scenarios. We es¬ 
tablish three new scenarios for comparison: three sensors and 30 targets, four sensors and 
40 targets, and five sensors and 50 targets. The number of sensors and targets maintain a 
constant ratio across scenarios; this allows us to observe the effect of preassigning sensors 
to targets on networks operating at their tracking capacity. Figure 4.4 displays the three 
sensor networks, including “best case” results from 60 simulations with and without parti¬ 
tioning. For six values of 7 , we obtain the “best case” schedule for 60 simulations of each 
of the three scenarios. 

Figure 4.5 and Figure 4.6 display the resulting delays, the number of observations, and the 
energy emitted. Again, the impact appears to be relatively stable until partitioning accounts 
for approximately 40 percent of sensor-to-target assignments. However, the relative impact 
of partitioning on the larger networks appear to be greater than that observed with three 
sensors. Worth noting are the average number of observations performed by the four- and 
five-sensor networks when the partitioning accounts for approximately 80 percent of target 
track assignments. Rather than continue the trend of requiring additional observations as 
more of the operating space is partitioned, these scenarios utilize less observations on aver¬ 
age than the scenarios where partitioning accounts for approximately 50 to 60 percent of the 
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target assignments. This is due to the signifieant inerease in delays for both networks when 
the partitioning method preassigns too many targets and overwhelms individual sensors. 

These preliminary results suggest that assigning a portion of targets to sensors based on 
their geographic location may provide an efficient heuristic approach for reducing the num¬ 
ber of real time sensor-to-target assignments that other, more sophisticated algorithms must 
perform. 
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Figure 4.4: Observations from optimal schedules, with and without geographic partition¬ 
ing, for 60 simulations of scenarios involving three, four, and five sensors. 
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Figure 4.5: Average delays, observations, and transmitted energy for scenarios involving 
three, four, and five sensors. 
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Figure 4.6: Delays, observations, and transmitted energy for each simulation of scenarios 
involving three, four, and five sensors. 
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4.3 Geographic Assignment with Handoff Penalties 

To account for resource costs resulting from the transfer of tracking responsibility between 
sensors, we utilize the formulation presented in Section 3.3. We consider the scenario 
where targets approach a single HVA from three directions. Figure 4.7 (a reproduction 
of Figure 3.5) represents the target density obtained by projecting Figure 3.2 onto two 
dimensions and discretizing the results into a 50 x 50 hexagonal grid. The three black lines 
represent the target paths for which penalties are assessed if a handoff between sensors 
occurs. During a target handoff, the sensor receiving responsibility for the target must 
perform extra work to establish the required track quality. Thus, the model only assesses a 
penalty on the sensor gaining the target. 



Figure 4.7: Two dimensional target density for one HVA with three directions of approach. 


To determine the potential effects of target handoffs on optimal geographic partitioning, we 
vary the handoff penalty a in 1 kW increments from 0 to 5 kW. Each instance contained 
approximately 30,000 decision variables, of which all but five were discrete, and 10,300 
constraints. Typical solution times were less than two seconds. Figure 4.8 and Table 4.1 
display the results of this analysis. 
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No Handoff Penalty 1 kW Handoff Penalty 



2 kW Handoff Penalty 


3 kW Handoff Penalty 




4 kW Handoff Penalty 5 kW Handoff Penalty 



Figure 4.8: Results of geographie assignment with various handoff penalties for the target 
density in Figure 4.7. 
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Handoff Penalty 
(a) 

kW 

Individual Sensor 

Sensor Requirements 

kW 

Total Network 

Sensor Requirements 

kW 

Number of 

Handoffs 

0 

3.97 

11.91 

2 

1 

4.97 

14.91 

2 

2 

5.53 

16.59 

1 

3 

5.67 

17.01 

1 

4 

5.79 

17.37 

0 

5 

5.79 

17.37 

0 


Table 4.1: Sensor energy demands and number of handoffs for geographic assignments 
with various handoff penalties for the target density in Figure 4.7. 


When there is no handoff penalty (a = OkW), the resulting partition matches that obtained 
from the continuous geographic partitioning formulation. In this scenario, even without 
handoff penalties, one of the target paths remains completely within the region assigned to 
Sensor 3 and no handoff occurs. As a increases to IkW, the other two target paths continue 
to transit between regions, but the partitioned regions adapt to the increased demands on 
the receiving sensor. Sensor 1. In order to minimize the maximum workload on any sensor. 
Sensor 1 shifts some of its region to other sensors, such as the lower right corner of the 
operating space, to equalize the increased burden across the sensor network. This transfer 
of regions between sensors costs the system additional resources since the range between 
the regions and their assigned sensors increases. Thus, while the total handoff penalty in 
this case is 2 kW, the total network radar resource requirement rises from 11.91 to 14.91 
kW, a 3 kW increase. 

When the handoff penalty exceeds approximately 1.5 kW, the costs of accepting two hand¬ 
off penalties exceeds the costs of reassigning regions among sensors and one of the target 
paths is completely absorbed by Sensor 1. As the handoff penalty increases to 2 and 3 kW, 
regions shift from Sensor 1 to equalize resource demands between sensors. 

Around 4 kW, the handoff penalty becomes so high that the sensor network does not accept 
any handoffs. In this scenario, each of the sensors maintains one of the target paths from 


51 




start to finish. Any additional increase in handoff penalty has no effect on the partitioned 
regions since no handoffs occur. 

While the partitioned regions adapt to any changes in handoff penalties, a reduction in the 
number of handoffs requires these penalties to be a large percentage of the network’s radar 
demands. The system removed a handoff only when the penalties exceeded approximately 
25 percent of total demands on the sensor network’s resources. To remove both handoffs, 
the total resource requirements for the sensor network increases by 45 percent. Thus, for 
the scenario considered, the amount of resources demanded to execute a handoff between 
sensors must be relatively high compared to other tracking requirements before it becomes 
advantageous to significantly alter the partitioned regions to limit handoffs. Analysis of 
additional scenarios may yield insights into the generality of these results. 
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CHAPTER 5: 
CONCLUSION 


The coordination of networked sensors is critical to obtaining superior functionality and 
battle space awareness. By synchronizing tracking efforts in a cluttered target environ¬ 
ment, individual platforms can redirect resources to unit-specific operations or other radar 
functions. This thesis develops several models to analyze target assignments in a multiple¬ 
sensor network. 

We develop an optimization model that uses full target information to provide a “best case” 
assessment of the ability of a given set of radar platforms to track a collection of targets. A 
modified version of this model determines the impact on these results if targets coordinate 
their maneuvers in an attempt to exceed the resource capacity of the sensor network. Re¬ 
sults suggest that synchronized maneuvers may produce situations whose sensor resource 
demands exceed network capabilities. This could result in reduced track quality leading to 
higher track loss rates. 

We also consider the more realistic scenario in which the planner’s knowledge is imperfect 
and describe approaches for allocating sensor assets to targets based on geographic loca¬ 
tion. The first method extends an existing two-dimensional geographic allocation approach 
to three dimensions. By adjusting how much of the operating space the model assigns to 
a specific sensor, we develop a heuristic approach for preassigning sensors to targets that 
has the potential to reduce the computational demands of more sophisticated assignment 
algorithms. In the scenarios examined, the heuristic can perform approximately 40 percent 
of sensor-to-target assignments with a limited impact on the “best case” assessment results. 
This could lead to significant reductions in real time computations and data transfer rates 
required by sensor-to-target assignment algorithms. 

Finally, we investigate the impact of transferring target tracking requirements between sen¬ 
sors by developing a model that permits the inclusion of handoff penalties. In order to do 
so, we modify the previously mentioned partitioning model by discretizing the operating 
space and defining target paths through this space. Any sensor that gains part of the path 
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after it begins is assessed a handoff penalty to account for the additional resources required 
when transferring a target between sensors. For the scenario considered, we find that signif¬ 
icant modifications to the sensors’ regional assignments only occur if the resource demands 
for target handoffs account for more than 25 percent of the total demand. 

5.1 Further Research Topics 

Several avenues of research exist to continue investigating this problem. A number of 
interesting extensions arise simply by relaxing the assumption that every observation results 
in the same information gain. 

Radar returns and retum-to-track correlation are not without error. An investigation into 
the effects of false negative and false positive observations along with the potential for mis- 
correlation could lead to a greater appreciation of actual resource requirements for a given 
scenario. A possible approach could include utilization of a dynamic programming model 
to analyze the potential impacts of choosing different waveforms or delaying observations. 

Environmental factors also have a significant impact on radar performance. Future analysis 
could explore this effect by adapting some of the models introduced in this thesis to allow 
for the emplacement of physical obstacles in the operating space, such as land or weather 
events. Additionally, target aspect can affect the RCS observed by different sensors. Ex¬ 
pansion of the models presented in this thesis to allow for unique target RCSs for each 
sensor is another possible research topic. 

In addition to external factors, design limitations also impact an MFR’s performance. MFR 
arrays are not capable of transmitting in all directions. As the angle between the direction 
of the radar beam and the array’s broadside increases, the emitted power becomes more 
diffuse and less effective at tracking. To overcome these limitations, the AN/SPY-1 MFR 
uses four arrays to provide 360 degree coverage. Future research could investigate how a 
sensor’s limited coverage angle affects the decision of how to schedule target observations. 

By allowing for varying levels of information gain, the possibility arises that the assignment 
of multiple sensors to a single target may result in superior track quality. This thesis does 
not investigate this possibility and the models presented provide a single sensor-to-target 
assignment. Development of a model that allows for multiple sensors-to-target assignment 
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during each time step would more closely represent the problem faced by networked MFRs. 

While relaxing our assumptions concerning information gain presents numerous challenges 
worth exploring, future research could also focus on the heuristic approach presented to 
allocate targets to sensors geographically. The proposed benefit of this method is to reduce 
the number of sensor-to-target assignment decisions that must be made in real time. Future 
analysis could determine the actual reduction in real time computations provided by this 
heuristic when implemented with existing assignment models, such as those presented in 
the literature review of Section 1.2. 

Future research could also compare the results to other geographic assignment models. 
It is possible that a simpler approach leads to similar or more effective sensor-to-target 
assignment. Additionally, the geometry of the sensor network might impact the effective¬ 
ness of the heuristic. If there are interior sensors, such as in a layered defensive position, 
only the exterior sensors have an advantage when observing distant targets. Limiting geo¬ 
graphic assignment to a subset of available sensors may result in superior performance of 
the heuristic. 

Several of the simulated scenarios result in periods of time during which targets overwhelm 
a single sensor’s tracking capability. To account for this, future research could focus on the 
development of a contingency algorithm that allows for the reassignment of targets that 
oversaturate a region. The algorithm would need to determine which targets become avail¬ 
able for other sensors to observe, when the targets become available, and how to transfer 
the responsibility for tracking between sensors. 

The challenges of protecting HVAs continue to grow with the rapid proliferation of aerial 
assets. Coordination of available defensive capabilities is critical to defeating this threat. 
This thesis provides several models for investigating the benefits of networked MFRs. 
As long as an effective air defense remains elusive, research of this topic will remain 
paramount. 
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Appendix: Proof of Duality and Optimality for 
Geographic Partitioning 


This appendix presents the derivation of the dual problem (Section 3.1.3) for the linear 
relaxation of the geographic assignment model presented in Section 3.1 and provides the 
proof of optimality for the partition presented as Equation (3.8). This appendix parallels 
work presented by Carlsson (2012) with modifications of the formulations to account for 
three dimensions and the specific parameters and nomenclature introduced in this thesis. 

Derivation of Dual to Geographic Assignment Model 

This derivation is a simple proof “sketch” for obtaining the dual formulation by discretizing 
the original problem. For a more rigorous mathematical proof, refer to Appendix A of 
(Carlsson, 2012). 

The geographic partition model considered is provided as Equations (1) and (2). 



( 1 ) 


( 2 ) 


This problem is reformulated as an infinite-dimensional integer program in Equations 3-6. 


min t (3) 

/i(-),-,/«(•) 

s.t. 
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Ii{x)aip{x)\\x - piW"^ dV 


V/ (4) 



D 


I4W = 1 

i=\ 

li{x)e {0,1} 


yxeD (5) 
Wi,x G D (6) 


where Ii{x) is 1 if ;c is assigned to faeility i and 0 otherwise. 

Without loss of generality, we can relax the equality constraint of Equation (5). Relaxation 
of the integer constraint (Equation (6)) results in the following infinite-dimensional linear 
program (Equations 7-10). 


min t (7) 

s.t. 

t> jjj cCip{x)\\x-pi\\^dV Vi (8) 

D 

n 

\/xeD (9) 

i=\ 

li{x)>0 \/i,xeD (\0) 


Discretization of the operating space into grid cells indexed by k admits the following 
formulation (Equations 11-14). 


mint (11) 

x.t 

s.t. 

t Vi (12) 

k 
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\fk (13) 


n 

^ 1 

i=l 

Xi,k>0 Vz,fc (14) 

where e is the volume of eaeh grid eell, Ci^k is the value of a,p (^) ||;c — evaluated at the 
eenter of grid eell k, and Xi^k represents the fraction of grid cell k assigned to sensor i. The 
dual of this formulation is presented as Equations 15-18. 


m&xYqk 

k 

(15) 

s.t. 


qk < 

Vz,zt (16) 

n 

i=l 

(17) 

ri,qk>0 

Vz,zt (18) 

By introducing the variable : 

— qk/£, we can rewrite the dual as Equations 19-22. 

maxe£4 

(19) 

k 


s.t. 


q'k < Ci,kri 

Vz,zt (20) 

i=\ 

(21) 

IV 

o 

(22) 


By reducing e so that the grid cells become infinitesimal, the dual formulation can be 
rewritten as Equations 23-26. 
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max jjj o{x)dV 

(23) 

D 


s.t. 


g{x) < ^iaip{x)\\x-piW"^ 

yx,i (24) 

VI 

(25) 

A,>0 

V/ (26) 


This is exactly the dual formulation presented in Section 3.1.3. 

Proof of Optimality of Geographic Partitioning Assignment 

This section provides the proof that Equation (3.9) (reproduced as Equation (27)) results in 
the optimal partitioning of the operating space for the model presented in Section 3.1.1. 

^,'=|jcgZ) I hcii\\x-pi\\^ <Xjaj'^x--y VyV*} (27) 

Theorem. Eet A* be the weight vector obtained by solving the dual problem (Equations 
23-26). Then setting 

^ 0 otherwise ’ 

with {R^j. ■ ■ = 9^(7),7, A*), is an optimal partition of D with respect to the infinite¬ 

dimensional integer programming formulation of the primal (Equations 3-6). 

Proof. Eor any nonzero objective value of the dual (Equation (23)), Equation (24) requires 
A/ > 0 for all i. Eor reasonable values (i.e., positive) of a,- and p(x), any feasible combina¬ 
tion of A,- > 0 for all i results in a positive objective value. Therefore the optimal objective 
value for the dual must be positive and A* > 0 for all i. It follows from complementary 
slackness that Equation (4) is binding for all i (t* = fff Ii{x)*aip{x)\\x — pi\f dV for all 

D 
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i). Thus,/// aip{x)\\x — pi\f dV = /// aip(^)||^ —jPill^dV for all/. Plugging A* and the 

R* R* 

o* (a:) into the dual yields 


JJJ a*{x)dV = IJJ minA*a,p(.r)||.r-p/||'^dy 

D ' b 

= ^ a/pWlk-P/f dy 

i=i 

= t^i III aiPWII^-Pifdy 

i — 1 D* 

n 

.r)||.r — Pill^dy (Because ^A* = 1) 

R* 

= jjj i;(x)a,p(x)\\x-pi\\'AV 

D 



which completes the proof. 
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